Quantum and Thermal Depinning of a String from a Linear Defect 
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I. INTRODUCTION 



Mikhail A. Skvortsov 
L. D. Landau Institute for Theoretical Physics, 117940 Moscow, Russia 
and Theoretische Physik, ETH-Honggerberg, CH-8093, Zurich, Switzerland 

The problem of a massive elastic string depinning from a linear defect under the action of a 
small driving force is considered. To exponential accuracy the decay rate is calculated with the 
help of the instanton method; then, fluctuations of the quasiclassical solution are taken into account 
to determine the preexponential factor. The decay rate exhibits a kind of first order transition 
from quantum tunneling to thermal activation with vanishing crossover region. The model may be 
applied to describe nucleation in 2-dimensional first order quantum phase transitions. 
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Macroscopic tunneling phenomena always attracted theorists' attention since it is the field where classical and 
t • quantum physics meet each other. It governs the low temperature decay of metastable states of physical systems with 
many degrees of freedom, separated from the neighboring states by high potential barriers; at high temperatures the 
metastable state decays via thermal activation. Several such systems were studied, including motion of dislocations 
across Peierls barriers , tunneling with strong friction jj| , or quantum breaking of an elastic string Q . 

In the present paper the problem of quantum and thermal depinning of a massive elastic string trapped in a linear 
defect and subject to a small driving force is considered. The interest in one-dimensional manifold tunneling has 
revived last years as it describes the creep of vortices trapped by columnar defects or plane twin boundaries ||. 
Unlike thermal activation, which for nondisspative Hamiltonian systems is independent of the dynamic properties of 
the string quantum tunneling is strongly dependent on the dynamics, which for the case of vortices in HTSC can 
be either dissipative or governed by a Hall force The tunneling problem of a massive string therefore cannot be 
^ • directly applied to the vortex creep but its study is important since it allows for an exact analytical solution. Massive 
string depinning from a linear defect describes the nucleation of a new phase in the vicinity of a first order quantum 
phase transition in 2 dimensions. The boundary of the two phases may be considered as a string, the difference in 
chemical potentials of the phases plays the role of a driving force. 

The problem is solved quasiclassically and exact solutions are obtained for the whole temperature range in the 
limit of small driving force. Then, fluctuations of the quasiclassical solutions are taken into account to determine the 
preexponentials. The special feature of the problem considered is the presence of a "first order" transition between 
the quantum tunneling and thermal activation regimes that should be contrasted with a smooth crossover between 
them in other related problems |j],|]J]| • 



II. MODEL AND SUMMARY 



Consider a massive elastic string on a plane pinned by a linear defect. The dynamics of the string is governed by 
the Lagrangian 



C[u(x, t)] 



dx 



(1) 



Here p is the mass density of the string, c is the speed of sound, and u is the displacement perpendicular to the string. 
The string is directed along the a;-axis. In the presence of a driving force, V(u) can be represented as 



V(u) = V (u) - Fu, 



(2) 



where Vq(u) is a potential of a single pinning line of size uq: Vq(0) — and Vb ( |u| 3> uq) = Eq. The driving force F 
renders the states localized near u — metastable and the string depins and leaves the defect in a finite time r —1 . 
The lifetime T _1 will be obtained analytically for the whole temperature range in the limit of small F. 
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At zero temperature the string has no energy to overcome the potential barrier that separates the states with u = 
and with u > Eq/F where the barrier disappears. In this case the decay rate is determined by macroscopic quantum 
tunneling. To study the tunneling process one has to rewrite the action in imaginary time (for temperature T): 



A[u(x,t)] 



h/2T 



H/2T 



(It I dx 



(3) 



and find the saddle point solution u(x, r) for this action that can be obtained from variation of ([}]). Such a solution 
of a Euclidean field theory is commonly referred to as an instanton ||. Then, to exponential accuracy the lifetime 
T" 1 oc exp(-A[u]/h). 

The action can be easily estimated as follows Let t be the tunneling time. Then the mass involved in the 

tunneling process can be estimated as M ~ per. We can describe the collective tunneling process as a tunneling of 
one particle of mass M through a potential barrier of height Eqct. The tunneling time can be obtained from the 
condition that the kinetic energy M(Eq/F) 2 /t 2 should be of the order of the potential energy Eqct. This gives the 
estimate r 2 ~ pEo j F 2 and the action 



PcE 2 
F 2 



(4) 



Such a quasiclassical treatment is valid provided that Aq/Ti 3> 1. This criterion can be expressed as 



< 1, 

pcu z 



(5) 



where a — F/F and F ~ E / u is the typical value of the force that keeps the string in the well. F is also of the 
order of the critical driving force F c at which the potential barrier disappears. The following consideration will be 
constrained to the limit of small a«l. 

At high temperatures when the periodicity h/T along the r-axis becomes smaller than the size ~ (pE ) 1/2 /F of the 
quantum instanton, the saddle point solution u(x) does not depend on r any more and the escape rate is determined 
by pure thermal activation. In this case T _1 oc exp(— A/T), where A is the energy of the optimal configuration u{x) 
that extremizes (|J). Note that it is independent of the dynamics since the thermal saddle point solution does not 
depend on imaginary time. The energy barrier A can be estimated in a way similar to the one used for the quantum 
case. Dimensional analysis yields 



A 



{pc 2 Elf/ 2 
F 



(6) 



The activation energy diverges at zero driving force. Such a behavior should be contrasted with the thermal kink- 
antikink nuclcation in the models with several degenerate minima (Sine Gordon, </> 4 theories) where in the 

absence of a driving force the activation barrier remains finite and is twice as large as the kink rest energy Ek ■ In 
the present case, the single minimum of Vq(u) is nondegenerate, resulting in large -F-dependent size of the instanton 
and divergent A. 

There exist a temperature T c when the probabilities of quantum and thermal instanton nucleations become equal, 



Jc ~% ~ (pE )i/ 2 - 



(7) 



Usually, this temperature indicates a crossover from quantum tunneling to thermal activation. This is the case for 
single particle tunneling j?j , tunnelling of a string between two degenerate minima when a small driving force lifting 
the degeneracy is applied [Qj^] , and quantum breaking of an elastic string . In the present problem the situation 
is different. It will be shown that for T < T c the string escapes due to quantum tunneling, temperature corrections 
to the action being exponentially small for q«1. For T > T c the lifetime is determined by pure thermal activation. 
Hence, at T — T c the system exhibits a kind of first-order phase transition from the quantum to the thermal regime 
and no crossover region in the conventional sense is obtained. Such a behavior is a consequence of the coexistence 
of two qualitatively different quantum and thermal instantons in the vicinity of T c . The total = (Tq + IV) -1 is 
equal to Tq 1 for T < T c and to T^ 1 for T > T c . On the other hand, in a usual situation there exist only one saddle 
point solution of (^) that describes quantum tunneling at low temperatures, pure activation at high temperatures, 
and crosses over between them at T ~ T r . 
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III. QUASICLASSICAL ANALYSIS 



In this section the saddle point solutions for the action (@) are obtained in the limit of small driving force. For 
further use I want to specify some quantities defined above by order of magnitude estimates. We will see below that 
in the limit of small F the quasiclassical results for the nucleation rate are independent of the certain features of 
Vq{u) at the scale Uq. The only relevant parameter is the depth Eq of the potential. The other characteristic of Vq(u) 
that determines instantons interaction and the strength of quantum fluctuations is the curvature at the origin Vq'(O). 
Instead of it I define the width of the pinning potential according to 



Un = 



2Eo_ 



The ratio of the driving force to the pinning force is the small parameter in the problem: 



Fup 



2F 2 



W(o) 



1/2 



In the following I will work in the units p = c = 1. One has to rescale x — > (pc 2 ) 
dimensional units. The classical equation of motion obtained from the action (||) is 

(dl + dl)u = V'{u). 



-V2, 



(8) 
(9) 

p _1 / 2 r to return to 
(10) 



The function u(x, r) is periodic in r with the period h/T and vanishes at |x| — > oo. Eq. ( |l0| ) defines a 2D nonlinear 
electrostatics problem, the charge density being dependent on the electrical potential u. 

The solution of Eq. ([To]) behaves essentially different for u <C Uq and for u ^> uq (see Eq. (0) for the definition of 
V(u)). Let us define a curve T on the (x, t) plane where u = uq. It separates the plane into inner and outer domains 
B and B. In the outer domain B, u < uq that corresponds to the string lying in the pinning well Vq(u). Inside the 
domain £>, u > uq and the string is under the barrier. 

The key point in the determination of u(x, r) from the nonlinear Eq. ( |Io|) is to obtain the solution in the outer 
domain B and to extract the boundary condition for the internal region of the instanton. The size of the instanton 
(R for the quantum instanton ([!(]), or L for the thermal instanton ([18]); see Fig. 1 for geometry clarifying) is of the 

1/2 

order of E ' / F as estimated above. In the external region B of the instanton, u rapidly falls off to zero at a distance 
a that can be estimated from Eq. (|i~(i| ) as a ~ [^"(O)] -1 / 2 ~ aR. 

Since a is much smaller than the boundary curvature radius (R for the quantum instanton, or oo for the thermal 
instanton), one can neglect the curvature in the Laplace operator in a layer of width a along the boundary T of the 
instanton and write simply V 2 ~ ~§^i where n is the outward normal to the boundary of the instanton. Then Eq. ( |l0| ) 
becomes the equation of motion of a classical particle in an inverted potential —V(u): 



(id 



The energy £ = \ ~ V( u ) is an integral of the motion. Its value should be obtained from the condition at infinity, 

that gives £ = 0. Applied to u = uq, the equation £ = gives the boundary condition for the interior problem: 



du 
dn 

The other condition is u|r = uq. Thus, we arrive at the following problem: 



V 2 u 



-F, 



du 
dn 



(12) 



(13) 
(14) 

(15) 



Eqs. ([l3]) and (|l4|) define a Dirichlet problem with a unique solution for a given contour T. The solution of Eq. ( |T3| ) 
with the boundary condition (^) (von Neumann problem) is unique as well. The problem (p^|)~(p^[) is therefore 
strongly overdetermined and cannot be solved for a general contour T. 
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There are only two regions for which the solution of Eq. ( |l3| ) satisfies both boundary conditions simultaneously. 
They are (i) the circle and (ii) the strip (see Fig. 1). 
(i) Circle of radius R. The solution is 



F 



u 



(R 2 ~r 2 ), 



where r 2 



t . The radius can be obtained from (12): 



R 2 



8E 

F 2 ■ 



(16) 



(17) 



This is the quantum instanton. 

(ii) Strip of width 2L parallel to the r-axis. Then 



F , 



u = u Q + ^ (I/ 2 



* 2 ) 



with 



L 2 = 



2£o 
F 2 ' 



(18) 



(19) 



This solution is time-independent and holds at large temperatures. 
The solutions are shown schematically in Fig. 2. 

To compute the action one has to know the solution of Eq. ([To]) everywhere. The contribution of the boundary 
layer to the action is not universal and depends on the features of the pinning potential Vq (u) but is of order a when 
compared to the contribution from the interior region. Thus, to leading order in a, the action for the circle is (here 
and below dimensional units are restored) 



An = 4tt 



pcEl 
F 2 ' 



whereas for the strip (temperature T) the result reads 

A T (T) 

with Ao the energy cost of the optimal tongue (|l8|), 

2 5/2 p l/2 cS 3/2 



A 
I . 

T ' 



Quantum and thermal actions become equal at 



3F 



T r . = 



2 1 / 2 hF 
3tt (pEo) 1 / 2 ' 



(20) 



(21) 



(22) 



(23) 



The strip ( |18| ) is a solution for all T. It will be shown in the next section that this solution becomes unstable below 
some temperature To = 0.9T C . The circle solution exists as long as its diameter 2R is less than the periodicity h/T in 
r. The quantum solution then ceases to exist at 



Ti = ^T c . 
8 

The dependence of the actions on temperature for the two solutions is illustrated in Fig. 3(a). 

The lifetime T -1 is determined by the solution with the minimal action. To exponential accuracy, 



(24) 



lnr- 



A Q /h for T<T C , 
A /T for T > T c . 



(25) 



This answer holds for all temperatures but not too close to T c . In a very narrow region around T c , where |1 — T/T c \ < 
a 2 — <C 1 (cf. (fib), both solutions contribute to T -1 according to the formula = (Fq + Tt)" 1 . Note 

P CU Q — 



T c /A 
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that the last equation is not a bridging formula but an exact expression up to exponentially small corrections. This 
is a consequence of the coexistence of two types of instantons close to T c . 

It is also instructive to study the temperature dependence of the quantum tunneling probability. It appears to 
be very small due to the exponential decay of u at large distances. Actually, the distance between the boundary of 
the instanton ( pf ) and its image is St = h/T — h/T±. Then according to Eq. (pp| ) the overlap of the exponentially 
decaying tails of the instanton and its image is as small as uq exp [-{V'^/py/^Sr}, resulting in 



An 



exp 



2E \ 1/2 /h h 



(26) 



In the vicinity of T\ , 



SAq(T) 
An 



exp 



T 



(27) 



IV. FLUCTUATIONS 



Let us now take fluctuations near the classical instantons into account and determine the preexponential factors in 
(p5|). The starting point is the expression for the lifetime |h1|| 



hr = 2T- 



^Ixsx Z 
ReZ' 



where Z is the partition function 



Z = y Du(x,r)exp 



A[u] 



(28) 



(29) 



The neighborhood of the saddle point solution u(x,t) of the equation 5A[u]/Su — gives the main contribution to 
the imaginary part of Z. I expand u near u according to 



u(x, t) = u(x, t) + 22 C a u a (x, t), 



(30) 



where u a (x,r) are the eigenfunctions of the operator 5 2 A[u}/Su 2 : 

{5 2 A[u]/8u 2 ) it u a {x,T) = X a u a (x,T). 



(31) 



The lowest eigenvalue Ao is negative, resulting in the imaginary contribution to the partition function. Ai = is 
also an eigenvalue of ( |3T[ ) related to the zero mode with respect to translation. For the thermal instanton (|T8|), Ai = 
is nondegenerate and u\{x) cx d x u corresponds to the zero mode along x-axis, while for the quantum instanton (|l(f), 
the zero level is twofold degenerate since the location of the instanton on the x, t plane depends on both coordinates. 

Applying standard methods [[n]Jl^|| for zero modes and the negative eigenvalues we arrive at the following expres- 
sions: The lifetime per unit length for thermal activation is 



r 

L 



T 
2-rrh 2 



dx 



<9tt x 
dx 



1/2 



det' (S 2 A[u]/5u 2 



det (S 2 A[u]/Su 2 ) l 



-1/2 



exp 



(32) 



The decay rate for quantum tunneling reads 

1/2 



r 
T 



l 

2nh 



dxdr ( — 
ox 



dxdr ( — — 

OT 



du N " 



1/2 



det" (S 2 A[u}/Su 2 ) u __ 



det (5* A[u]/Su 2 ) u=0 

Here det' and det" denote the determinants with the zero eigenvalues omitted. 



-1/2 



exp 



A[i 



(33) 
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A. High-temperature instanton 



For the high temperature r-independent solution ( [1 8|) (I use the units h = p = c = 1), 

S 2 A 



^=-d 2 -d 2 x + U(x); (34) 



^) = ^». (35) 



The eigenvalues of the operator (34) are given by 

A m „ = e m + (2irTn) 2 , (36) 
where n is an integer and {e m } is the spectrum of the one dimensional Schrodinger equation 

[-dl + U(x)]tp m (x) = e m ip m (x). (37) 



The potential U(x) is illustrated schematically in Fig. 4. L oc oT x is the semiwidth of the instanton strip (p"8|), the 
width of the minimum of U(x) near \x\ — L, a ~ La, is independent of a; U(\x\ ^ L) — uo 2 = Vq'(O) is the oscillator 
frequency for the pinning potential Vq(u). The form of the potential near |x| = L depends on the detailed features 
of Vq(u). For a — ► 0, L — ► oo and the potential well C/(x) becomes deep in the sense that as many as -a -1 discrete 
levels exist in this potential. Excited states can be obtained within the quasiclassical approximation. It is possible 
also to obtain the low lying levels in the limit of small a. 

Let us define 7 as a function of the continuous parameter e as 

7 (e) =d x \n^{L), (38) 

where VK 2 -) is the solution of Eq. (|37]) with energy e that vanishes at x — ► 00. Once the function 7(e) is known, one can 
obtain the discrete spectrum in the following way. For |x| < L, the solutions of Eq. ( p7| ) are given by ipn( x ) — cos(^/ex) 
and ip e {x) = sin(Y / ex) for even and odd states respectively. The spectrum then can be obtained from the continuity 
of the logarithmic derivative of ipe( x ) &t x = L: 

d x ln^ e (L)=7(e). (39) 

The function 7(e) is not universal and depends on the shape of the minimum of U(x), that is on the details of Vq(u). 
It appears however that 7(0) is independent of the pinning potential. We may find 7(0) by using our knowledge that 
the first excited level has zero energy. For \x\ < L its wave function linearly depends on x yielding 

7 (0) = l/L. (40) 



The low lying states can be obtained from Eq. (|39|) with the right hand side set equal to 7(0) thus neglecting the 
difference between 7(e) and 7(0) for small e. Then one immediately obtains the implicit equations for the spectrum: 

ViX tan( VeX) = -1 (41) 

for even states and 

y/lLcot^L) = I (42) 

for odd states. 

The approximation 7(e) ~ 7(0) is valid for small enough energies e <C e* ~ 7(0)/7'(0). One can estimate 7'(0) as 
^{lo 2 )/lo 2 ~ l/w, yielding e* ~ aT 1 / L 2 . On the other hand, the energy of the n-th excited level can be estimated as 
e„ ~ (jm/ \2L)) 2 for large n. Therefore, Eqs. (41) and ( [42]) correctly describe as many as cr 1 ! 2 low lying levels. 

The negative eigenvalue £q can be obtained from Eq. Mil) and is equal to 



eo = -j^, (43) 



where fi = 1.19968 is the solution of the equation ^ttanh/i = 1. 



6 



The negative eigenvalue determines the temperature To where Aoi changes the sign leading to the instability of the 
thermal instanton ( |l8| ) with respect to a small perturbation of its boundary. Inserting ( f43| ) into Eq. (|36|), I obtain 

in dimensional units, 

hl 0.9007;. (45) 



" 2 3 / 2 7T ( P E y/ 2 

Now we are in a position to calculate the ratio of determinants in (j32|). To do it, one has to compute the products 
over n 0], resulting in 



det' (S 2 A[u]/6u 



2 ) - 

I U—U 



det (S 2 A[u]/5u 2 ) u=0 
with the quantum correction to the barrier height, 



-1/2 



7 e- Al/T x(^), (4<i) 



2sin^ A Vr 



A - 1 V* ,V2 , /'°° rf g ^(q) 1/2 r47 X 
m=2 y 

Here (5(g) is the phase shift of the states in the continuous spectrum, e(q) — oj 2 + q 2 , and M labels the last discrete 
state. The function \ is dehned according to 

InX ( f) = - £ ln(l - e-^/ T) _ f ln(1 _ e -e ( ,)V 2/T) _ (48) 



Tn 7 ^— i 7 _ 2tt dq 

*> ' m=2 J-oo I 

At large q, S(q) can be calculated perturbatively JljJ, 

8{q)=-±- J[U{x)-u 2 ]dx = ^, (49) 

leading to a logarithmic divergence of the integral over continuous spectrum in ( |47| ) that should be cut off at some 
wave vector A: 

u 2 L A 

Ai = -— In-. 50 

47T W 

The ultraviolet cutoff A is a new parameter that should be introduced for the problem described by the Lagrangian 
(|l|) to be well defined. The value of A remains undefined in the model in question. Physically, it is determined by the 
length scale at which the elastic approximation of the bending energy in (|l]) breaks down. Then A may be estimated 
as the inverse core radius of the linear manifold considered. 

x{T/T c ) is a complicated function of temperature. For T close to To, only the lowest level with positive energy e 2 
is excited. Its energy can be obtained from Eq. ([51): €2 = (2.79839/L) 2 , resulting in 

( 14.7\ 

X (x~l)~l + expf —J. (51) 

For intermediate temperatures, To -C T <C To/a ~ ui, many levels contribute to the sum in (|4§|). In this limit, 
e m ~ (nm/(2L)) 2 and 

X (x) = cxp . (52) 

For T 3> To, many levels are excited and the quasiclassical approach can be applied. The sum over the states in 
Eq. (Eq) can be expressed via an integral over phase space: 



7 



lnx 



-21n(l 



) + J_ [ dxdp L(l _ e -V^/T ) _ ln(1 _ e -y/u&H?IT) 



-21n(l -e 



2tt 

-T /Tn , ^£ 



dp 



ln(l - e -V^?/ T ) - ln(l - e~ p / T ) 



(53) 



The first term in fl53| ) accounts for the two omitted states in Eq. (E§|). For T <C T Q /a, Eq. ( |53| ) reduces to Eq. (|52 
For T > T /a, one finds 



(54) 



Substituting (111) and @ into @, we obtain 



T /A 
I =B(T)exp -- 



(55) 



A = Aq + Ai and B(T) in conventional units are given by the following expressions: 



A 



2 5/ 2(0 l/2 c£; 3/2 



3h 



3F 



8irpciiQ 



In 



Au Q 



PC 

Eo 



2\ 1/2 



(56) 



B(T) 



p l/4 E 3/i T 3/2 , r 



2 l/4 3 l/2^1/2 tf c l/2pl/2 ginI Tb X I Tq 



(57) 



B. Quantum instanton 



Fluctuations of the quantum instanton (16) are governed by the following Hamiltonian 



^ = -a»-fl2 + CT(r); 



where /7(r) is cylindrically symmetrical, 



U(r) = 



d 2 V{u{r)) 
du 2 



(58) 



(59) 



The potential U (r) looks very similar to the one-dimensional potential U(x > 0) illustrated in Fig. 4. 
The ratio of determinants in (|3^) can be expressed in terms of the spectrum as 



dot" (5 2 A[u]/Su 2 



det (S 2 A[u]/Su 2 ) u=0 



-1/2 



exp 



|-^E" ln l e ™l + ^E ln l £ ™l f ' 



(60) 



where the first sum is taken over the states in the potential (|5^) and the second in the potential U (r) = lu 2 . 

We can estimate (^0|) within the quasiclassical approach by taking advantage of the large number of bound states 
in the potential ( |59| ) for small a. According to the quasiclassical quantization rule, each state occupies the volume 
(27r) 2 in phase space and (BOf) reduces to 



eexp 



1 1 



d 2 rd 2 p 



In y/U(r) + p 2 - In ^/u 2 + p 2 



2 (2tt) 2 

The integral over the momentum is logarithmically divergent, leading to 

'R 2 uo 2 , A 



e exp 
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pdphi 



uj 2 + p 2 



eexp 



ln- 

LU 



(61) 



(62) 



Inserting (p21) in (p3), and using (Km rewritten in conventional units, I obtain the final result 
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The results ( p6|) and 
provides the condition 



r 
I 



■ exp 



A = Air 



F 2 



In 



1/2' 



(63) 



(64) 



are valid provided that the fluctuation contribution is less than the classical one. This 



h 



2irpcuQ 



In 



Au 



PC 

E 



1/2' 



< 1. 



(65) 



This is a much stronger condition than that of (||) . When (^35|) fails then one has to take unharmonic fluctuations into 
account that will renormalize the numerical coefficient in front of (|64|), whereas the estimate (^) holds as long as the 
condition (||) is fulfilled. 



V. CONCLUSION 



The problem of a massive string depinning from a linear defect has been solved in the limit of small driving force 
(a <C 1) allowing for an exact analytical solution in the whole temperature range. It has been shown that there 
exists a temperature region T < T < T\ where the action (^) has two different saddle point solutions: ( |l6| ) and 
(|l8|). The actions for these solutions become equal at T = T c . For T < T c , the lifetime is determined by quantum 
tunneling (Eq. ([63])) with exponentially small temperature corrections (p7|). For T > T c , the lifetime is given by a 
pure activation expression (|55|). At T — T c the system jumps from quantum to classical behavior. The dependence 
of lnT -1 vs. T at constant F is shown in Fig. 3(a). The possibility of such a first order transition between quantum 
and thermal behavior was first discovered in Ref. H] , 

It is worth emphasizing that the first order transition is obtained for small a <C 1. As a increases, To (a) and 
Ti(a) come closer to each other and become equal at some a c ~ 1, that is, there exists a tricritical point where the 
first order transition disappears. For a close but below a c , the temperature corrections to the quantum tunneling 
rate become large and the transition at T c (a ) should be considered as a transition from thermally assisted quantum 
tunneling to pure activation. 

Besides the transition at constant F, one can consider the transition at constant T by changing the applied force 
which may be easier to achive experimentally. lnT -1 as a function of F is shown in Fig. 3(b). For F < F c = 
■^jij-^Eo) 1 / 2 , the metastable state decays via thermal activation with lnT -1 cx F^ 1 , whereas for F > F c , the string 
escapes via quantum tunneling with lnT -1 cx F~ 2 . 

The problem in question is connected with the nucleation phenomenon in 2-dimensional first order quantum phase 
transitions. Consider a first order phase transition controlled by some parameter x. The quantum nature of the 
transition implies that it exists even at T = 0. One possible example is the liquid-solid transition in 2D helium films. 
Let x > correspond to phase I, and x < to phase II. The phase boundary is a string that has mass density p 
and elasticity pc 2 — a, where a is the boundary line tension. The difference in chemical potentials of the phases 
is equivalent to the driving force F = fii — flu oc x. Suppose that phase II is prepared at x > 0. Such a state is 
metastable and is destroyed by the nucleation of phase I that can take place near the edges of the 2D system or in 
the bulk. Nucleation near the edges is exactly the problem of a massive string depinning from a linear defect. In 
the present case, the edge of the system plays the role of the pinning well with the depth Eq = a. It can be shown 
that nucleation always occurs near the edges since the bulk tunneling action and activation barrier are larger than 
the corresponding values for the edge: Ab u ik/A e d ge — 37r/2 5 / 2 and A-bulk/ Aedge = 16\/2/15. 

Finally, the results ( |56| ) and ( |64] ) for the string on a plane can be easily generalized to the string in (1 + d) dimensions. 
The instantons (|l6| ) and (|l^) remain the same since any transverse displacement will cost extra bending energy. The 
preexponentials will be different; they will be multiplied by an additional factor from transverse fluctuations. The 
main effect is the appearance of a coefficient d in front of the logarithm in formulae ([56]) and (^J) . 



ACKNOWLEDGMENTS 



I am grateful to M.V.Feigel'man for his support and fruitful discussions, to G. Blatter, V.B.Geshkenbein, L.S.Levitov 
and D.M.Kagan for helpful comments. I express my gratitude to the Institute for Theoretical Physics, ETH-Ziirich, 



9 



Switzerland, for hospitality. This work was supported by the joint grant M6M300 from the ISF and the Russian 
Government, by the RFBR grant 95-02-05720 and by the Swiss National Foundation. 



[1] B. V. Petukhov and V. L. Pokrovsky, Zh. Eksp. Teor. Fiz. 63, 634 (1972) [Sov. Phys. - JETP 36, 336 (1973)]. 
[2] B. I. Ivlev and V. I. Mel'nikov, Phys. Rev. B 36, 6889 (1987). 

[3] A. I. Larkin and Yu. N. Ovchinnikov, Zh. Eksp. Teor. Fiz. 86, 719 (1984) [Sov. Phys. - JETP 59, 420 (1984)]. 
[4] L. S. Levitov, A. V. Shytov, and A. Yu. Yakovets, Phys. Rev. Lett. 75, 370 (1995). 
[5] For a review, see G. Blatter et al, Rev. Mod. Phys. 66, 1125 (1994). 

[6] As an example of classical nucleation with dissipation, see F. Marchesoni, Phys. Rev. Lett. 73, 2394 (1994) and references 
therein. 

[7] I. Affleck, Phys. Rev. Lett. 46, 338 (1981); U.Weiss, Quantum Dissipative Systems, World Scientific (1993). 

[8] R. Rajaraman, Solitons and Instantons, North-Holland Publishing Company, Amsterdam/New York/Oxford (1982); S. 
Coleman, Aspects of Symmetry, Cambridge University Press (1985). 

[9] M. I. Dyakonov, Fiz. Tver. Tela 29, 2587 (1987). 
[10] J. S. Langer, Ann. Phys. 41, 108 (1967). 
[11] I. Zittartz and J. S. Langer, Phys. Rev. 148, 741 (1966). 
[12] C. Callan and S. Coleman, Phys. Rev. D 16, 1762 (1977). 

[13] L. D. Landau and E. M. Lifshitz, Quantum Mechanics, Course in Theoretical Physics Vol. 3, Pergamon, London/Paris 
(1958). 

[14] I. M. Lifshitz and Yu. Kagan, Zh. Eksp. Teor. Fiz. 62, 385 (1972) [Sov. Phys. - JETP 35, 206 (1972)]. 



10 



FIG. 1. The geometry of the quantum (a) and thermal (b) instantons. 



FIG. 2. The form of the thermal instanton. u(x) is given by (18) for \x\ < L and exponentially decays for \x\ > L. The 
figure can also be considered as the radial dependence u(r) of the quantum instanton (16). 

FIG. 3. (a) The dependence of the actions of the quantum (Qu) and thermal (Th) instantons on temperature at constant 
F. The decay rate exhibits a transition from quantum tunneling to thermal activation at T — T c . (b) Decay rate vs. F at 
constant T. F = F c separates the quantum and thermal regimes. 

FIG. 4. Potential energy for the problem (37) linearized in the vicinity of the thermal instanton. Also shown are the wave 
functions of the ground (dashed line) and the first excited (dotted line) states. 



11 



1 


'u 


-L 


L x 



Fig. 2 





U(x)' 


or 


\ 

s 

— ► 






l\ 


X 



Fig. 4 



